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Abstract 

On the basis of the confinement mechanism earlier proposed by author the electric 
formfactor of 7r°-meson is nonperturbatively calculated. The latter is then applied to 
describing electromagnetic decays ir° — ► 27 and 7r° — ► e + e~7 which entails estimates 
for parameters of the confining SU(3)-gluonic field in 7r°-meson. The corresponding 
estimates of the gluon concentrations, electric and magnetic colour field strengths 
are also adduced for the mentioned field. 
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1 Introduction and preliminary remarks 



In Refs. [1-3] for the Dirac- Yang-Mills system derived from QCD-Lagrangian 
there was found and explored an unique family of compatible nonperturba- 
tive solutions which could pretend to decsribing confinement of two quarks. 
Applications of the family to description of both the heavy quarkonia spectra 
[4] and a number of properties of pions and kaons [5] showed that the confine- 
ment mechanism is qualitatively the same for both light mesons and heavy 
quarkonia. 

Two main physical reasons for linear confinement in the mechanism under 
discussion are the following ones. The first one is that gluon exchange be- 
tween quarks is realized with the propagator different from the photon one 
and existence and form of such a propagator is direct consequence of the 
unique confining nonperturbative solutions of the Yang-Mills equations [2,3]. 
The second reason is that, owing to the structure of mentioned propagator, 



Preprint submitted to Elsevier 



1 February 2008 



quarks mainly emit and interchange the soft gluons so the gluon condensate (a 
classical gluon field) between quarks basically consists of soft gluons (for more 
details see Refs. [2,3]) but, because of that any gluon also emits gluons (still 
softer), the corresponding gluon concentrations rapidly become huge and form 
the linear confining magnetic colour field of enormous strengths which leads 
to confinement of quarks. This is by virtue of the fact that just magnetic part 
of the mentioned propagator is responsible for larger portion of gluon con- 
centrations at large distances since the magnetic part has stronger infrared 
singularities than the electric one. Under the circumstances physically nonlin- 
earity of the Yang- Mills equations effectively vanishes so the latter possess the 
unique nonperturbative confining solutions of the Abelian-like form (with the 
values in Cartan subalgebra of SU(3)-Lie algebra) [2,3] that describe the gluon 
condensate under consideration. Moreover, since the overwhelming majority 
of gluons is soft they cannot leave hadron (meson) until some gluon obtains 
additional energy (due to an external reason) to rush out. So we deal with 
confinement of gluons as well. 

The approach under discussion equips us with the explicit wave functions for 
every two quarks (meson). The wave functions are parametrized by a set of real 
constants aj,bj,Bj describing the mentioned nonperturbative confining gluon 
field (the gluon condensate) and they are nonperturbative modulo square inte- 
grable solutions of the Dirac equation in the above confining SU(3)-field and 
also depend on /i , the reduced mass of the current masses of quarks forming 
meson. It is clear that under the given approach just constants aj,bj, Bj, fj, 
determine all properties of any meson and they should be extracted from ex- 
perimental data. 

Such a program has been to a certain extent advanced in Refs. [4,5]. The 
aim of the present paper is to continue obtaining estimates for aj,bj,Bj for 
concrete mesons starting from experimental data on spectroscopy of one or 
another meson. We here consider 7r°-meson and its electromagnetic decays 
7r° — > 27 and 7T° — > e + e~7 whose widths practically add up to the full width 
of 7r°-meson [6]. 

Of course, when conducting our considerations we shall rely on the standard 
quark model (SQM) based on SU(3)-flavor symmetry (see, e. g., Ref. [6]) so 
in accordance with SQM ir° ={uu — dd)/y/2 is equiprobable superposition of 
two quarkonia, consequently, we shall have two sets of parameters a,j,bj,Bj. 

Further we shall deal with the metric of the flat Minkowski spacetime M that 
we write down (using the ordinary set of local spherical coordinates r, 1?, (p for 
the spatial part) in the form 

ds 2 = g^dx* ® dx v = dt 2 - dr 2 - r 2 (d$ 2 + sin 2 tidy 2 ) . (1) 

Besides, we have \8\ = \det{g lxl> )\ = (r 2 sini9) 2 and 0<r<oo, 0<-#<7r, 
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< ip < 2ir. 



Throughout the paper we employ the Heaviside-Lorentz system of units with 
h — c—1, unless explicitly stated otherwise, so the gauge coupling constant g 
and the strong coupling constant a s are connected by relation g 2 /(47r) = a s . 
Further we shall denote L 2 (F) the set of the modulo square integrable complex 
functions on any manifold F furnished with an integration measure, then 
LQiF) will be the n-fold direct product of L 2 (F) endowed with the obvious 
scalar product while f and * stand, respectively, for Hermitian and complex 
conjugation. Our choice of Dirac 7-matrices conforms to the so-called standard 
representation and is the same as in Ref. [5]. At last ® means tensorial product 
of matrices and I3 is the unit 3x3 matrix. 

When calculating we apply the relations 1 GeV -1 « 0.1973269679 fm, 1 s _1 
0.658211915 x 10" 24 GeV , 1 V/m « 0.2309956375 x 10~ 23 GeV 2 , 1 T « 
0.6925075988 x 10~ 15 GeV 2 . 

Finally, for the necessary estimates we shall employ the T 00 -component (volu- 
metric energy density ) of the energy- momentum tensor for a SU (3)- Yang-Mills 
field which should be written in the chosen system of units in the form 



2 Specification of main relations 

As was mentioned above, our considerations shall be based on the unique fam- 
ily of compatible nonperturbative solutions for the Dirac- Yang-Mills system 
(derived from QCD-Lagrangian) studied in details in Refs. [1-3]. Referring for 
more details to those references, let us briefly decribe and specify only the 
relations necessary to us in the present Letter. 

One part of the mentioned family is presented by the unique nonperturbative 
confining solution of the Yang-Mills equations for A = A^dx^ = A a ^\ a dx^ (A a 
are the known Gell-Mann matrices, /i — t, r, ■&, <p, a — 1, ...,8) and looks as 
follows 



Tp, = -F a m F a v , g«P + \F^ i^tf W • 



(2) 




+ A 1 , An = 




2 • 



Am = 




ai + a 2 



r 



(A-i + A 2 ) , 



A\ip 




bir + Bi , A 2ip 




b 2 r + B 2 , 



3 



A 



3ip 



jLa 8 

V3 * 



-{b 1 + b 2 )r-{B 1 + B 2 



(3) 



with the real constants aj, Aj, bj, Bj parametrizing the family. As has been re- 
peatedly explained in Refs. [2-5], parameters A^ 2 of solution (3) are inessential 
for physics in question and we can consider A\= A 2 = 0. Obviously we have 
Y%=i Ajt = Z)f=i = which reflects the fact that for any matrix T from 
SU(3)-Lie algebra we have TrT = 0. 

Another part of the family is given by the unique nonperturbative modulo 
square integrable solutions of the Dirac equation in the confining SU(3)-field of 
(3) \& = (^i, ty 2 , ^3) with the four-dimensional Dirac spinors fyj representing 
the jth colour component of the meson, so \& may describe relative motion 
(relativistic bound states) of two quarks in mesons and the mentioned Dirac 
equation looks as follows 




#1 














H 2 




(::)■ 


H 2 ^ 2 








H 3 J 







(4) 



where Hamiltonian Hj is 



Hi 



1 



/i - i^dr -i'j 2 - [ d# 
r 



2 / r sin v 



-91° (fAjt 



+ r 



^—A 



9^ + - sin^7 1 7 3 + - cos-#7 2 7 3 



(5) 



with the gauge coupling constant g while Hq is a mass parameter and one can 
consider it to be the reduced mass which is equal, e. g., for quarkonia, to half 
the current mass of quarks forming a quarkonium. 

Then the unique nonperturbative modulo square integrable solutions of (4) 
are (with Pauli matrix <Ji) 



e -iwjt r -l 



F^r)*^) 
F j2 {r)a^,ip) 



,3 = 1,2,3 



(6) 



with the 2D eigenspinor $~ 



.72 



of the euclidean Dirac operator D 011 



the unit sphere S 2 , while the coordinate r stands for the distance between 
quarks. The explicit form of $j is not needed here and can be found in Refs. 
[3,7]. For the purpose of the present Letter we shall adduce the necessary 
spinors below. Spinors form an orthonormal basis in Ll(§ 2 ). We can call 
the quantity ujj relative energy of jth colour component of meson (while ipj 
is wave function of a stationary state for jth colour component) but we can 
see that if we want to interpret (4) as equation for eigenvalues of the relative 
motion energy, i. e., to rewrite it in the form Hip = uoip with = (ijji, ip 2 , ijj 3 ) 
then we should put u = ujj for any j so that Hjipj = ojjipj = uipj. Under this 
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situation, if a meson is composed of quarks qi j2 with different flavours then 
the energy spectrum of the meson will be given by e = m qi + m q2 + u with 
the current quark masses m Qk ( rest energies) of the corresponding quarks. On 
the other hand for determination of oj 3 the following quadratic equation can 
be obtained [1-3] 

\g 2 a 2 +(n 3 +a 3 ) 2 }uj 2 -2(\ 3 -gB 3 )g 2 a 3 b 3 u j +[(X j -gB j ) 2 -(n j +a j ) 2 ]g% 2 j - f j 2 (n j +a j ) 2 = 0, 

(7) 

that yields 

(jj j — UJ j [71 j , Ijj \ j ' ^ — 

A 3 g 2 a 3 b 3 ± (n 3 + a^y/jn] + 2n 3 a 3 + A 2 )fi 2 + g 2 b 2 (n 2 + . _ 

n 2 + 2n 3 a 3 + A 2 ,J - 1,2,3, 

(8) 

where a 3 = -(ai + a 2 ), h = -(h + b 2 ), B 3 = -(B 1 + B 2 ), A 3 = X 3 - gB 3 , 
a 3 = \J A 2 — g 2 a 2 , Uj = 0, 1, 2, while \j = ±(l 3 + 1) are the eigenvalues of 
euclidean Dirac operator V on unit sphere with l 3 = 0, 1,2, .... It should be 
noted that in previous papers [1-5] we used the ansatz (6) with the factor e tuJjt 
instead of e~ tuJjt but then the Dirac equation (4) would look as —id t ^ = 
and in equation (7) the second summand would have plus sign while the first 
summand in numerator of (8) would have minus sign. We here return to the 
conventional form of writing Dirac equation and this slightly modifies the 
equations (7)-(8). In line with the above we should have oo = uj\ = u 2 = u 3 in 
energy spectrum e = m qi +m q2 +u for any meson and this at once imposes two 
conditions on parameters a 3) b 3) B 3 when choosing some experimental value for 
e at the given current quark masses m qi ,m q2 . 

Within the given Letter we need only the radial parts of (6) at n 3 = (the 
ground state) that are 

f 31 = c,iy-< (i - ^r) > Pi = 9b 3 + ft, 

F j2 = K -,()/-< ■ (l + 9 -^j , Qj = & - u s (9) 

with (3j = yj nl — uj 2 + g 2 b 2 while C 3 is determined from the normalization con- 
dition jo°°(l-Fji| 2 + \F 32 \ 2 )dr = |. Consequently, we shall gain that G -^2( r3 ) 
at any i G R and, as a result, the solutions of (6) may describe relativistic 
bound states (mesons) with the energy (mass) spectrum e. 

It is useful to specify the nonrelativistic limit (when c — > oo) for spectrum 
(8). For that one should replace g — > g/y/hc, a 3 — > aj/Vhc, b 3 — > bjVhc, 
Bj — > Bj/y/hc and, expanding (8) in z = 1/c, we shall get 

uJ 3 (n 3 J 3 ,X 3 ) = 



5 



±/x c 



it 



2h\n j + |Aj 



+ 



\jg 2 djbj 
K n i + l A il 



=F A*o 



g 3 B j a 2 j f(n j , \j) 
h^rij + lXjiy 



where f(rij, Xj) 



AXjUj^ + X^ + ^^ + ^ + Xf). 



z+0(z 2 ) 
(10) 



As is seen from (10), at c — > oo the contribution of linear magnetic colour 
field (parameters bj,Bj) to spectrum really vanishes and spectrum in essence 
becomes purely Coulomb one (modulo the rest energy). Also it is clear that 



when rij — > oo, 0;^ — > ±^JJc^~+~g 2 b 



We may seemingly use (8) with various combinations of signes (±) before 
second summand in numerators of (8) but, due to (10), it is reasonable to 
take all signs equal to plus which is our choice within the Letter. Besides, 
as is not complicated to see, radial parts in nonrelativistic limit have the 
behaviour of form Fji, Fj2 ~ r lj+1 , which allows one to call quantum number 
lj angular momentum for jth colour component though angular momentum 
is not conserved in the field (3) [1,3]. So for meson under consideration we 
should put all lj — 0. 

Finally it should be noted that spectrum (8) is degenerated owing to degen- 
eracy of eigenvalues for the euclidean Dirac operator V on the unit sphere 
§ 2 . Namely, each eigenvlalue of V X = ±(l + 1), I — 0, 1, 2..., has multiplicity 
2(1 + 1) so we has 2(1 + 1) eigenspinors orthogonal to each other. Ad referen- 
dum we need eigenspinors corresponding to A = ±1 (I = 0) so here is their 
explicit form 



A = 



A 




i<fi/2 



c 



, or $ = — 



e 2 



2 V -e 



or * = 

2 V e l 2 



-i<p/2 



-iip/2 



with the coefficient C = l/y/2n (for more details see Refs. [3,7]). 



(11) 



3 Electric and magnetic formfactors, anomalous magnetic moment 



As has been mentioned in Section 1, we shall analyse electromagnetic decays 
of 7r°-meson so let us consider what quantities characterizing electromagnetic 
properties of a meson we could construct within our approach. 

Within the present Letter we shall use relations (8) at rij = = lj so energy 
(mass) of meson under consideration is given by /i — 2m q + uj with oj = 
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ujj(0, 0, Xj) for any j = 1, 2, 3 whereas 

g 2 a i b 1 ai/xo g 2 a 2 b 2 , a 2 /i 2 a 3 &3 , «3/Uo /1oA 

^ = — l \- T1TT = — a \- T1TT = — a \- T1TT = H~ 2m <? (12) 

Ai |Ai| A 2 |A 2 | A 3 |A 3 | ^ 9 

and, as a consequence, the corresponding meson wave functions of (6) are 
represented by (9), (11). 



3.1 Choice of quark masses and gauge coupling constant 



It is evident for employing the above relations we have to assign some val- 
ues to quark masses and gauge coupling constant g. In accordance with Ref. 
[6], at present the current quark masses necessary to us are restricted to in- 
tervals 1.5 MeV < m u < 3 MeV, 3.0 MeV < m d < 7 MeV, so we take 
m u = (1.5 + 3)/2 MeV = 2.25 MeV, m d = (3 + 7)/2 MeV = 5 MeV. Under 
the circumstances, the reduced mass of (5) will respectively take values 
m u /2,m<i/2. As to gauge coupling constant g = y / 4~7raj, it should be noted 
that recently some attempts have been made to generalize the standard for- 
mula for a s = a s (Q 2 ) = 127r/[(33 — 2nf) In (Q 2 /A 2 )] (rif is number of quark 
flavours) holding true at the momentum transfer \[Q 1 ~ — > oo to the whole 
interval < y/Q 2 < oo. We shall employ one such a generalization [8] used in 
Refs. [9]. It looks as follows (x = ^ in GeV) 

a[X) ~ (33-2^)1^+/^ 1 J 

with 

A(o:) = l+^ (1 + 3:)( 1 3 2 3 ' 2n/) ln^-l) 1 + 0.6x 13 j , f 2 (x) = m 2 (l+2.8s 2 )- 2 , 

wherefrom one can conclude that a s — > 7r = 3.1415... when x — > 0, i. e., 
g _> 2tt = 6.2831.... We used (13) at m = 1 GeV, A = 0.234 GeV, n f = 3, 
x = m n o = 134.976 MeV to obtain g = 6.101482007 necessary for our further 
computations at the mass scale of 7r°-meson. 



3. 2 Electric formf actor 



For each meson with the wave function ^/ = (\I/j) of (6) we can define electro- 
magnetic current J* = <g> I 3 )V = (^ t ^, ^(a <g> 7 3 )#) = (p, J), a = 7%. 
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Electric formfactor f{K) is the Fourier transform of p 

f(K) = f &$e- iKT d 3 x = E / ^je~ iKr d 3 x = E/jW 

J n=l J 7=1 



-iKr 



E I {\F n \ 2 + \F j2 \ 2 )&]&j?—2-cPx, d 3 x = r 2 witidrd&dip (14) 

with the momentum transfer K. At rij = = as is easily seen, for 
any spinor of (11) we have = l/(47r), so the integrand in (14) does 

not depend on ip and we can consider vector K to be directed along z- 
axis. Then Kr = Kr cos-# and with the help of (9) and relations (see Ref. 

[10]): f^r^e-^dr = T(a)p~ a , Re a,p > 0, ^r a ~ l e^ r (^(^r)) dr = 

r(a)(^+p 2 )-/ 2 f ^S a ^SK), Re a > -1, Rep > |Imif|, r(a + 
V cos (a arctan (K/p)) ) 

1) = aiT (a), Jo e- iKrcos *smi)dtf = 2 sin (Kr)/(Kr), we shall obtain 
f(K , ^ f( K\ f ( 2 ^) 2aj+1 sin [2a, arctan (if/^))] 

wherefrom it is clear that /(if) is a function of if 2 , as should be, and we can 
determine the root-mean-square radius of meson in the form 



<r >- 



y 2af + 3a 3 - + 1 



It is clear, we can directly calculate < r > in accordance with the standard 

quantum mechanics rules as < r >= yj J r 2 ^^d 3 x = ]J I r 2 ^^jd 3 x and 

the result will be the same as in (16). So we should not call < r > of (16) 
the charge radius of meson - it is just the radius of meson determined by the 
wave functions of (6) (at n 3 - — = lj) with respect to strong interaction. Now 
we should notice the expression (15) to depend on 3-vector K. To rewrite it 
in the form holding true for any 4- vector Q, let us remind that according to 
general considerations (see, e.g., Ref. [11]) the relation (15) corresponds to the 
so-called Breit frame where Q 2 = —K 2 [when fixing metric by (1)] so it is not 
complicated to rewrite (15) for arbitrary Q in the form 



3 a/g/j ^+i S i n [2 a , arctan (J\Q 2 \/(2f3 j ))] 

f(Q 2 ) = E W 2 ) = E ^ ^ ^ V 1 J (17) 

3=1 j=l 0a 3 

which passes on to (15) in the Breit frame. 
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3.3 Magnetic moment, magnetic form] "actor, anomalous magnetic moment 



We can define the volumetric magnetic moment density by m = q(r x J)/2 = 
Q_[{yJz — z Jy)i + {zJx — xJ z )i + (xJ y — yJ x )k]/2 with the meson charge q and 
J = ^(a <S> h)^- Using (6) we have in the explicit form 

3 $ 3 &0~20~i<&j 

J x = ^{F* 1 F j2 + F* 2 Fji)^-, J y = ^(F* l F j2 - F* 2 Fji)- 3 



^ = E(W-W^^ as) 

with Pauli matrices 01,2,3 ■ Magnetic moment of meson is M = J v md 3 x, where 
V is volume of meson (the ball of radius < r >). Then at rij = lj = 0, as is 
seen from (9), (11), F* x = Fji,F* 2 = —Fj 2 , < lv-020i ( £j = for any spinor of 
(11) which entails J x = J y = 0, i.e., m z = while J v m X;y d 3 x = because 
of turning the integral over ip to zero, which is easy to check. As a result, 
magnetic moments of mesons with the wave functions of (6) (at lj = 0) are 
equal to zero, as should be according to experimental data [6]. 

We can, however, define magnetic formfactor F(K) of meson if noticing that 



J \Jm x + m 2 + m 2 d 3 x ^ 
v 



( J m x d 3 x) 2 + ( J m y d 3 x) 2 + ( J m z d 3 x) 2 — M — 0. 
v v v 



Under the circumstances we should define F(K) as the inverse Fourier trans- 



form of J ml + m 2 + m 2 



m 2 x + m 2 y + m 2 z = — j F(K)e iKr d 3 K 



2/z(2tt) s 

with meson mass \i so that 



F(K) = ^ J yfi% + m 2 y + m 2 e- tKr d 3 x = fi J r \ J z sin i? | e~ iKr d 3 x 



^ f pWW^-e^x (19) 

for any spinor of (11) so the integrand in (19) does not again depend on (p. Then 
Kr = Kr cos-# and using the relation Jq e - lKrcosi) s i n 3 fid-d = 4[sin (Kr)/(Kr) 3 - 
cos (Kr)/(Kr) 2 }, we shall obtain 



oo cj 

F(K) = AfM frJ2C]P,Q 

n .7=1 



,2a, -2 PjT I , 9 2 b]\ [sin (Kr) cos (Kr) 







(Kr) 3 (Krf 



dr 
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2/i (2f3 j ) 2a ^ +1 PjQA 1 -^-) 



3 <Xj (1 - f )2 P 2 + (1 + 9 _hy Q 2 K 2 {K 2 + 4/3 2 )Q , 

- cos [2a, arctan (K/(2/3j))] \ = 



sin [(2a j - 1) arctan (K/(2(5 j ))\ 



[2aj -l)K{K 2 + 4ff)- 1 / 2 



3 

2 u2 



-a 3 {2a 3 + 1) a,-(4aj + 12a 2 + 11a, + 3) r2 + Q 



3 ft (i - f) 2 ^ 2 + (i + f) 2 ^ 2 L 6 ^ 120 ^. 

(20) 

where the fact was used that by virtue of the normalization condition for wave 
functions we have C 2 [P 2 (l-gbj/ fof+Q^l+gbj/ fy) 2 ] = (2/3,) 2 ^ +1 /[3r(2a,+ 
1)]. It is clear from (20) that F(K) is a function of K 2 and we can rewrite 
(20) for arbitrary 4-vector Q by the same manner as was done for electric 
formfactor f(K) in (17). 

Now we can define the anomalous magnetic moment density for meson by the 
relation 

so anomalous magnetic moment is 

M a = j m a d 3 x « J m a d 3 x = ^-F(0) (21) 

V R3 

with 

F(0) ~ ^ V 2 ^ + 1 Pjgj ' (1 - 

3ft 6/3, (1 - S*t)2 P 2 + (1 + ^ )2 q2 > 

as follows from (20). It is clear that for neutral mesons M a = due to q — 0. It 
should be noted that a possibility of existence of anomalous magnetic moments 
for mesons is permanently discussed in literature (see, e.g., Refs. [12] and 
references therein) though there is no experimental evidence in this direction 
[6]. 



4 Estimates for parameters of SU(3)-gluonic field in 7r°-meson 



4-1 Basic equations 



The question now is how to apply the obtained formfactors to electromagnetic 
decays n° — > 27 and n° — > e + e~7 to estimate parameters of SU(3)-gluonic 
field in 7r°-meson. The corresponding diagrams for decays under consideration 
are presented in Fig. 1. Let us at first consider two-photon decay (Fig. la). 



10 



Actually kinematic analysis based on Lorentz- and gauge invariances gives rise 
to the following expression for width F 1 of the given decay (see, e.g., Refs. [13]) 



ri = l^Lgiy (22) 



with electromagnetic coupling constant a em =l/137.0359895 and 7r°-meson 
mass pi = 134.976 MeV while the information about strong interaction of 
quarks in 7r°-meson is encoded in a decay constant g nil . Making replacement 
g njJ = y/2/ (4ir 2 f P ) we can reduce (22) to the more conventional form (see, e. 
g., Refs. [14]) 

r '=Mf (23) 

so the present-day experimental value I\ ps 7.741 eV conforms to fp 
130 MeV [6]. Further parametrization in the form f P = /z<3> entails $ ~ 
0.9673759912. We can now notice that the only invariant which <£> might 
depend on is Q 2 = fi 2 , i. e. we should find such a function f(Q 2 ) for that 
f(Q 2 = fi 2 ) ~ 0.9673759912. It is obvious from physical point of view that / 
should be connected with electromagnetic properties of 7r°-meson. As we have 
seen above in Section 3, there are at least two suitable functions for this aim 
- electric and magnetic formfactors. But, as was mentioned, there exist no 
experimental consequences related to magnetic formfactor at present whereas 
electric one to some extent determines, e. g., an effective size of meson in the 
form < r > of (16). It is reasonable, therefore, to take just electric formfactor 
of (17) for the sought function. Then, denoting the quantities fi/(2(3j) = Xj, 
we obtain the following equation for parameters of the confining SU(3)-gluonic 
field in 7r°-meson 

HQ* =M 2 )=E/i(Q 2 = M 2 )=E ^V ^^'^^ » 0.9673759912. 

j=l j=l Ua j X j \ L X j) 

(24) 

One more equation can be gained from consideration of one-photon decay (Fig. 
lb) with width T 2 ~ 0.0939 eV. The latter in crucial way depends on the value 
of electric formfactor of 7r°-meson at Q 2 = (2m e ) 2 , where m e = 0.510998918 
MeV is electron mass (for more details see Ref. [6] and book by Pilkuhn in 
[13]). The present-day value for f(Q 2 = Am 2 e ) is between 0.9999936936 and 
1.000001835 [6]. So, when designating m e / ' j3j = i/j, we obtain the mentioned 
equation in the form 

m - 4<, - 1 - - 1 • sin( ^;r fe) - -o. 

(25) 

At last, equations (12) and (16) should also be added to (24), (25) and the sys- 
tem obtained in such a way should be solved compatibly if taking < r >=0.672 
fm [6]. 
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Fig. 1. Two- and one-photon electromagnetic decays of tt -meson 
4-2 Numerical results 

The results of numerical compatible solving of equations (12), (16), (24) and 
(25) are adduced in Tables 1-3. 

Table 1 



Gauge coupling constant, mass parameter hq and parameters of the confining SU(3)- 
gluonic field for 7r°-meson 



Particle 


9 


Mo (MeV) 


dl 


0.2 


h (GeV) 


b 2 (GeV) 


#1 


B 2 


7T° UU 


6.10148 


1.125 


-0.0434737 


-0.00680835 


0.0848234 


0.0433136 


0.01 


-0.150 


7r —dd 


6.10148 


2.50 


-0.0606679 


0.0251427 


0.0956303 


0.0648174 


0.1250 


-0.2450 



Table 2 

Theoretical and experimental 7r°-meson mass and radius 



Particle 


Theoret. fi (MeV) 


Expcrim. n (McVJ 


Theoret. < r > (fm) 


Expcrim. < r > (fm) 


7T° UU 


a = 2m u + ojj(0, 0, -1) = 134.976 


134.976 


0.674927 


0.672 


TT°—dd 


fi = 2m d +wj(0,0,-l) = 134.976 


134.976 


0.675676 


0.672 



Table 3 

Theoretical and experimental 7r°-meson electric formfactor values 



Particle 


Theoret. /(Q 2 = p 2 ) 


Expcrim. f(Q 2 = fj 2 ) 


Theoret. /(Q 2 = 4m 2 ) 


Experim. /(Q 2 = 4m 2 ) 


7T° UU 


0.991172 


0.967376 


0.999999464 


0.9999936936-1.000001835 


TT°—dd 


0.996480 


0.967376 


0.999999762 


0.9999936936-1.000001835 



5 Estimates of gluon concentrations, electric and magnetic colour 
field strengths 

Now let us remind that, according to Refs. [3,5], one can confront the field (3) 
with Too-component (volumetric energy density of the SU(3)-gluonic field) of 
the energy-momentum tensor (2) so that 

_ _ rr _ E 2 + H 2 _ 1 (aj + a^ + al bj + b x b 2 + b\\ _ A B 

-tOO = J-tt — ~ — 77 7 1 2 • 2 q = ~\ + o • 2 q 

2 2 y r 4 r z sin v J r 4 r A sin v 

(26) 

with electric E and magnetic H colour field strengths and with real A > 0, 
B > 0. 

To estimate the gluon concentrations we can employ (26) and, taking the 
quantity u — T, the full decay width of a meson, for the characteristic fre- 
quency of gluons we obtain the sought characteristic concentration n in the 
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form 



_ Too_ 



(27) 



so we can rewrite (26) in the form T 0Q = Tqq 111 + Tqq 1 conforming to the con- 
tributions from the Coulomb and linear parts of the solution (3). This entails 
the corresponding split of n from (27) as n = n cou \ + n\ in . 

The parameters of Table 1 were employed when computing and for simplicity 
we put sin-$ = 1 in (26). Also there was used the following present-day full 
decay width of 7r°-meson: r = 1/r with the life time r = 8.4 x 10~ 17 s, whereas 
the Bohr radius a = 0.529177249 • 10 5 fm [6]. 

Table 4 contains the numerical results for n cou i, nn n , n, E, H for the meson 
under discussion. 

Table 4 



Gluon concentrations, electric and magnetic colour field strengths in ir -meson 



7T° — uu: 


r =< r >= 0.674927 fm 






















r 

(fm) 


n coul 






"lin 






n 




E 

(V/m) 






H 

(T) 




0.1r 


0.893407 X 


10 57 




0.232356 


X 


10 56 


0.916643 X 


10 57 


0.174836 


X 


10 24 


0.476625 X 


10 15 


ro 


0.893407 x 


10 53 




0.232356 


X 


10 54 


0.321696 X 


10 54 


0.174836 


X 


1Q 22 


0.476625 X 


10 14 


1.0 


0.185386 X 


10 53 




0.105844 


X 


10 54 


0.124383 X 


10 54 


0.796426 


X 


1() 21 


0.321687 X 


10 14 


10ro 


0.893407 X 


1Q 49 




0.232356 


X 


10 52 


0.233249 x 


10 52 


0.174836 


X 


10 20 


0.476625 X 


10 13 


a 


0.236413 X 


10 34 




0.377976 


X 


10 44 


0.377976 X 


10 44 


0.284409 


X 


1() 12 


0.607900 X 


10° 


vr — Id: 


r =< r ; 


>= 0.6756' 


7 6 fm 






















r 

(fm) 


n coul 
(m- 3 ) 






"lin 
(m- 3 ) 






n 

(m- 3 ) 




E 

(V/m) 






H 

(T) 




0.1r 


0.111063 X 


10 S8 




0.355535 


X 


10 56 


0.114618 X 


1Q 58 


0.194936 


X 


10 24 


0.589577 x 


10 15 


ro 


0.111063 X 


10 54 




0.355535 


X 


10 54 


0.466597 X 


10 54 


0.194936 


X 


10 22 


0.589577 x 


10 14 


1.0 


0.231485 X 


10 53 




0.162315 


X 


10 54 


0.185464 x 


10 54 


0.889956 


X 


10 21 


0.398363 X 


10 14 


l()r 


0.111063 X 


10 s " 




0.355535 


X 


10 52 


0.356645 X 


10 52 


0.194936 


X 


10 20 


0.589577 x 


10 13 


a O 


0.295201 x 


10 34 




0.579638 


X 


10 44 


0.579638 X 


10 44 


0.317809 


X 


1() 12 


0.752797 x 


10° 
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6 Discussion and concluding remarks 

6. 1 Discussion 

As is seen from Table 4, at the characteristic scales of 7r°-meson the gluon con- 
centrations are huge and the corresponding fields (electric and magnetic colour 
ones) can be considered to be the classical ones with enormous strenghts. The 
part n coul of gluon concentration n connected with the Coulomb electric colour 
field is decreasing faster than rii in , the part of n related to the linear mag- 
netic colour field, and at large distances ri\\ n becomes dominant. It should be 
emphasized that in fact the gluon concentrations are much greater than the 
estimates given in Table 4 because the latter are the estimates for maximal 
possible gluon frequencies, i.e. for maximal possible gluon impulses (under the 
concrete situation of 7r°- meson). The given picture is in concordance with the 
one obtained in Refs. [4,5]. As a result, the confinement mechanism developed 
in Refs. [1-3] is also confirmed by the considerations of the present Letter. 

It should be noted, however, that our results are of a preliminary character 
which is readily apparent, for example, from that the current quark masses (as 
well as the gauge coupling constant g) used in computation are known only 
within the certain limits and we can expect similar limits for the magnitudes 
discussed in the Letter so it is neccesary further specification of the parameters 
for the confining SU(3)-gluonic field in 7r°-meson which can be obtained, for 
instance, by calculating weak formfactors in processes K + — > 7r° + l + + v\ ( 
I — e, /i) with the help of wave functions discussed above. We hope to continue 
analysing the given problems elsewhere. 

At last it should be noted that the relation (23) is often connected with the 
so-called chiral (axial) anomaly and PCAC model (see for more details, e.g., 
Refs. [14]). Under the latter approach, however, the value of fp in (23) is 
all the same taken from experimental data while we express it through more 
fundamental parameters connected with the unique exact solution (3) of Yang- 
Mills equations describing confinement. But, e. g., in the case of decay rj — > 27 
the PCAC model is unlikely to be applicable and our approach works since it 
is based on the unique family of compatible nonperturbative solutions for the 
Dirac- Yang-Mills system derived from QCD-Lagrangian and, as a result, the 
approach is itself nonperturbative and may be employed for any meson. 

6.2 Concluding remarks 

Finally we should note the following. Our approach actually repeats the one of 
quantum electrodynamics (QED) at large distances: nonperturbative effects 
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of quark interaction are described by the modulo square integrable solutions 
of Dirac equation in the classical confining SU(3)-gluonic field (condensate of 
huge number of gluons) and for this aim we use the unique compatible non- 
perturbative solutions of the Dirac- Yang-Mills system derived directly from 
QCD-Lagrangian which is in perfect analogy with QED where one uses the 
unique modulo square integrable solutions of Dirac (Schrodinger) equation 
in the Coulomb field (condensate of huge number of photons), i. e., one em- 
ploys the unique compatible nonperturbative solutions of the Maxwell-Dirac 
(Schrodinger) system directly derived from QED-Lagrangian to describe, for 
example, the hydrogen atom spectrum. But the same analogy with QED 
prompts to us that further study should perhaps include radiative correc- 
tions. This should, however, be conditioned by experiment. At the moment, 
as far as is known to us, little is known about such effects. At any rate, as 
seems to us, our approach treats nonperturbative meson physics in a sound 
way. Another matter that QCD-Lagrangian probably allows one to develop 
some other approaches (not based on compatible solutions of the Dirac- Yang- 
Mills system) to confinement but our one appears to be rather natural since 
it is practically identical to the standard approach of quantum mechanics and 
QED to description of bound states in external (electromagnetic) fields. In 
other words, this approach should have been one of the very first approaches 
as soon as the QCD-Lagrangian was written. Historically, however, this way 
was rejected due to incomprehensible reasons without any detailed study. 
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